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We present updated measurements of time-dependent CP asymmetries in fully-reconstructed 
neutral B decays to several CP eigenstates containing a charmonium meson. The measurements 
use a data sample of (383 ± 4) x 10 6 T(4S) -» BB decays collected with the BABAR detector 
at the PEP-II B factory. We determine sin2/3 = 0.714 ± 0.032 (stat) ± 0.018 (syst) and |A| = 
0.952 ± 0.022 (stat) ± 0.017 (syst). 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



The standard model (SM) of electroweak interactions 
describes CP violation as a consequence of an irreducible 
phase in the three-family Cabibbo-Kobayashi-Maskawa 
(CKM) quark- mixing matrix In the CKM frame- 
work, neutral B decays to CP eigenstates containing a 
charmonium and a meson through tree-diagram 

dominated processes provide a direct measurement of 
sin2/3 Q, where the angle j3 is defined in terms of the 
CKM matrix elements V tJ as ^g[-(V cd V; b )/(V td V t * b )}. 

We report updated measurements, based on a sample 
of (383 ± 4) x 10 6 T(45) -> BB decays, of sin2/3 and 
of the parameter |A|. Here A = (q/p)(A/A) [§], q and 
p are complex constants that relate the P-meson flavor 
eigenstates to the mass eigenstates, and A/ A is the ratio 
of amplitudes of the decay of a B° or B° to the final 
state under study. We reconstruct B° decays to the fi- 
nal states J/ipK°, J/ipK°, ^(2S)K°, XdK°, VcK° and 
J/ip K*° 0] . Since our previously published result Q , we 
have added 157 x 10 6 BB decays and applied improved 
event reconstruction algorithms to the entire dataset. We 
have also developed a new r] c K® event selection based on 
the Dalitz plot structure of the ?y c — ► K®K + tt~ decay, 
and have performed a more detailed study of the CP 
properties of the background events, which results in re- 
duced systematic errors. We now include the J/ipK® and 
J/ipK*° modes in the sample to measure |A|, and we re- 
port individual measurements of sin2/3 and |A| for each 
of the CP decay modes used in the analysis. Finally, we 
present separate results for the J/ipK^-K+n^ +7r°7r°) 
and J/iPK°(K° + K°) modes. 

We identify (tag) the initial flavor of the reconstructed 
B candidate, B lcc , using information from the other B 
meson, Stag, in the event. The decay rate g+ (g~) for a 
neutral B meson decaying to a CP eigenstate accompa- 
nied by a B° (B°) tag can be expressed as 



S±(At) 



22mA 



,-\At\/Tg 

4r R o 



sin(AmrfAt) 



(1 =F Aw) ± (1- 2w) x 
cos(Am t iA<) 



1 + A 



(1) 



where At = t ICC — f tag is the difference between the proper 
decay times of the reconstructed and tag B mesons, rgo is 
the neutral B lifetime and Am c i is the mass difference of 
the B meson mass eigenstates determined from B° — B° 
oscillations 0] ■ We assume that the corresponding decay- 
width difference ATd is zero. The average mistag prob- 
ability w describes the effect of incorrect tags, and Aw 



is the difference between the mistag probabilities for B° 
and B° . The sine term in Eq. [T] results from the inter- 
ference between direct decay and decay after B° — B° 
oscillation. A non-zero cosine term arises from the in- 
terference between decay amplitudes with different weak 
and strong phases (direct CP violation) or from CP viola- 
tion in B°—B° mixing. In the SM, CP violation in mixing 
and direct CP violation in b — > ccs decays are both neg- 
ligible Q. Under these assumptions, A = r//e _24/3 , where 
rjf = ±1 is the CP eigenvalue of the final state /. Thus, 
the time-dependent CP-violating asymmetry is 



g + {At)+g^(At) 
-(1 - 2w)r) } sin2/?sin(Am d At). 



(2) 



The BABAR detector is described in detail elsewhere 
We select a sample of neutral B mesons (Bcp) de- 
caying to the r\j — —1 final states J/ipK®, ip(2S)K®, 
XciK® and rj c K®, and to the r/f = +1 final state J/ipK®. 
We reconstruct K® — > 7r + 7r~, except in J/ipK®, where 
we also include K® — > 7r°7r°. The charmonium mesons 
are reconstructed in the decays J/ip — > e + e~, n + fj.~; 
ip{2S) — > e + e~, fi + fjT , J/ipTT + n~ ; \ci — * J/^l an d 
?y c -> K°K + tt~. We also reconstruct the J/ipK*°(K*° -> 
K®ir°) final state, which can be CP-even or CP-odd due 
to the presence of even (P=0, 2) and odd (L=l ) orbital 
angular momentum contributions. Ignoring the angular 
information in J/ipK*° results in a dilution of the mea- 
sured CP asymmetry by a factor 1 1 — 2R± | , where R± is 
the fraction of the L=l contribution. In Ref. [j| we have 
measured R±_ = 0.233 ± 0.010 (stat) ± 0.005 (syst), which 
gives an effective r\ f = 0.504 ± 0.033 for / = J/ijjK* , 
after acceptance corrections. 

In addition to the CP modes described above, we use a 
sample of B° mesons (B& av ) decaying to the flavor eigen- 
states DW"ft+ (h+ = 7r+,p+,a+) and J/ipK*° (K*° -» 
K + tt~) to calibrate the flavor-tagging performance and 
At resolution. We also perform studies to measure ap- 
parent CP violation arising from CP-conserving processes 
using a control sample of B + mesons decaying to the fi- 
nal states J/ipK^+, ^{2S)K+ , Xc\K + , and r) c K+ . The 
event selection and candidate reconstruction remain un- 
changed from those described in Rcfs. SEE!, with the 
exception of modes containing r\ c mesons. In Ref. Q we 
reconstructed the B° — > ri c K® and P ± — > r\ c K ± modes 
using the r\ c — > K®K + ir~ decay, with the requirement 
2.91 < m K o K+7T - < 3.05GeV/c 2 . We now exploit the 
fact that the r\ c decays predominantly through a Kix res- 
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onance at around 1430MeV/c 2 and a K®K resonance 
close to threshold, and require that either m K o^- or 
m K+TT- be in the mass-range [1.26, 1.63] GeV/c 2 , or that 
m K+K o G [1.0,1.4] GeV/c 2 . 

We calculate the time interval At between the two B 
decays from the measured separation Az between the de- 
cay vertices of B lec and -Btag along the collision (z) axis 
[101 ] . The z position of the B rec vertex is determined from 
the charged daughter tracks. The B tag decay vertex is de- 
termined by fitting tracks not belonging to the B lec can- 
didate to a common vertex, while employing constraints 
from the beamspot location and the B lec momentum [Io| . 
Events are accepted if the calculated At uncertainty is 
less than 2.5 ps and |At| is less than 20 ps. The fraction 
of all events satisfying these requirements is 95%. 

The algorithm used to determine the flavor of the B tag 
at its decay to be either B° or B° is described in detail 
in Ref. [H . In brief, we define six mutually exclusive tag- 
ging categories in order of decreasing tag purity: Lepton, 
Kaon I, Kaon II, Kaon-Pion, Pion and Other. The figurc- 
of-merit for tagging is the effective tagging efficiency 
Q = J2i £ i(^ ~ 2wi) 2 , where Si is the tagging efficiency 
of tagging category i. We measure Q = (30.5 ± 0.3) %, 
consistent with the results in Ref. [f|. 

We determine the composition of our final sam- 
ple using the beam-energy substituted mass tties = 

V(^bcam) 2 - (P* B ) 2 > Where f beam and P*B are the beam 

energy and B momentum in the e + e center-of-mass 
(CM) frame. For the J/ipK® mode we instead use the 
difference AE between the candidate CM energy and 
£beam- The composition of our final sample is shown 
in Fig. [1] We use events with m E s > 5.2GeV/c 2 
{\AE\ < 80McV for J/ipK°) to determine the proper- 
ties of the background contributions. We define a signal 
region 5.27 < m ES < 5.29GeV/c 2 (\AE\ < lOMcV for 
J/i/jK®), which contains 12677 CP candidate events that 
satisfy the tagging and vertexing requirements (see Ta- 
ble HI). For all modes except r] c K° and J/ipK°, we use 
simulated events to estimate the fractions of events that 
peak in the jjjes signal region due to cross-feed from other 
decay modes (peaking background). For the r] c K® mode, 
the cross-feed fraction is determined from a fit to the 
wkk-k and Toes distributions in data. For the J/ip K® de- 
cay mode, the sample composition, effective 77/ , and AE 
distribution of the individual background sources are de- 
termined either from simulation (for B — > J/ipX) or from 
the mg+e- sidebands in data (for non- J/ip background). 

We determine sin2/3 and |A| from a simultaneous max- 
imum likelihood fit to the At distribution of the tagged 
Bcp and Bflav samples. The At distributions of the Bqp 
sample are modeled by Eq. [TJ Those of the -Bfl av sample 
evolve according to Eq. [T]with A = 0. The observed am- 
plitudes for the CP asymmetry in the Bcp sample and 
for flavor oscillation in the i?fl av sample are reduced by 
the same factor, 1 — 2w, due to flavor mistags. The At dis- 
tributions for the signal are convolved with a resolution 
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m ES (GeV/c 2 ) m ES (GeV/c 2 ) 

FIG. 1: Distributions for Bcp and _Bfl av candidates satisfying 
the tagging and vertexing requirements: a) tties for the final 
states J/tpK° s , ip{2S)K° s , Xd#°, and r) c K° s , b) AE for the 
final state J/ip K° , c) m ES for J/ipK*°(K*° -> K%ir°), and d) 
mES for the Bfl av sample. In each plot, the shaded region is 
the estimated background contribution. 

function common to both the Bfl av and Bqp samples, 
modeled by the sum of three Gaussian functions [13] ■ The 
combinatorial background is incorporated with an empir- 
ical description of its At spectra, containing prompt and 
non-prompt lifetime components convolved with a reso- 
lution function [10( distinct from that of the signal. The 
peaking background is assigned the same At distribution 
as the signal but with no CP violation, with the same At 
resolution function. 

In addition to sin2/3 and |A|, there are 68 free param- 
eters in the CP fit. For the signal, these are the pa- 
rameters of the At resolution (7), the average mistag 
fractions w and the differences Aw between B° and B° 
mistag fractions for each tagging category (12), and the 
difference between B° and B° reconstruction and tag- 
ging efficiencies (7). The background is described by 
mistag fractions (24), parameters of the At resolution (3) 
and Bflav time dependence (3), and parameters for the 
CP background (8), including the apparent CP asymme- 
try of non-peaking events in each tagging category. Fi- 
nally, we allow for the possibility of direct CP violation 
in the \ciK® background to J/ipK*° (1), and in the main 
backgrounds to the J/ipK® mode, coming from J/ipK®, 
J/ipK*°, and the remaining J/ip background (3 param- 
eters). The effective |A| of the non- J/i[> background is 
fixed from a fit to the J/tp -candidate sidebands in J/ipK®. 
We fix t b o = 1.530 ps and Am d = 0.507 ps" 1 Q. The 
determination of the mistag fractions and At resolution 
function parameters for the signal is dominated by the 
.Bflav sample, about 10 times more abundant than the 



TABLE I: Number of events iVtag and signal purity P in 
the signal region after tagging and vertexing requirements, 
and results of fitting for CP asymmetries in the Bcp sample 
and various subsamples. In addition, fit results for the Bfl av 
and B + control samples demonstrate that no artificial CP 
asymmetry is found where we expect no CP violation (sin2/3 = 
0, |A| = 1). Errors are statistical only. 



MOO 



Sample 


N tag 


P(%) 


sin2/3 


|A| 


Full CP sample 


12677 


75 


0.714 ±0.032 


0.952 ± 0.022 


J/lPKl (7T+7T-) 


4459 


96 


0.702 ± 0.042 


0.976 ± 0.030 


JH>K°g (W ) 


1086 


88 


0.617 ±0.103 


0.812 ±0.058 


■4>{2S)K° S 


687 


83 


0.947 ±0.112 


0.867 ± 0.079 




313 


89 


0.759 ±0.170 


0.804 ± 0.102 


VcK° s 


328 


69 


0.778 ±0.195 


0.948 ± 0.141 


JM>kI 


4748 


55 


0.734 ± 0.074 


1.061 ± 0.063 


J/ipK*° 


1056 


66 


0.477 ± 0.271 


0.954 ± 0.083 



J/iPK 


10275 


76 


0.697 ± 0.035 


0.966 ± 0.025 


J/i>K° s 


5547 


94 


0.686 ± 0.039 


0.950 ± 0.027 


Vf = -l 


6873 


92 


0.711 ±0.036 


0.935 ± 0.024 


1999-2002 data 
2003-2004 data 
2005-2006 data 


3084 
4850 
4725 


79 
77 
74 


0.735 ± 0.063 
0.728 ± 0.052 
0.681 ±0.054 


0.987 ± 0.045 
0.940 ± 0.035 
0.940 ± 0.037 



Lepton 
Kaon I 
Kaon II 
Kaon-Pion 
Pion 
Other 



1349 80 

1843 76 

2948 72 

2321 73 

2551 76 

1665 73 



0.728 
0.689 
0.751 
0.654 
0.671 
0.705 



± 0.066 
± 0.063 
± 0.071 
± 0.112 
±0.167 
± 0.504 



0.901 
0.986 
0.880 
0.999 
0.927 
1.506 



± 0.043 
± 0.046 
± 0.044 
± 0.075 
± 0.104 
± 0.483 



Bfl av sample 


123893 


85 


0.018 ±0.010 


0.995 ± 0.007 


B + sample 


29598 


94 


0.012 ±0.017 


1.010 ±0.012 



CP sample. 

The fit to the Bcp and £?ft av samples yields sin2/3 = 
0.714 ± 0.032 and |A| = 0.952 ± 0.022, where the errors 
are statistical only. The correlation between these two 
parameters is —1.5%. We also perform a separate fit in 
which we allow different sin2/3 and |A| values for each 
charmonium decay mode, a fit to the J/ipK® (ir + ir~ ± 
7r°7r°) mode, and a fit to the J/ipK° {K° + K°) sample. 
We split the data sample by run period and by tagging 
category. We perform the CP measurements on control 
samples with no expected CP asymmetry. The results of 
these fits are summarized in Table HI The difference in 
the 77 c -ftT2sin2/3 value with respect to our previous pub- 
lication [5j is partly due to the slightly different recon- 
struction algorithms and partly to the different selection; 
the two measurements are consistent when the system- 
atic error is taken into account. 

Figure [2] shows the At distributions and asymmetries 
in yields between events with B>° tags and B° tags for 
the 77/ = —1 and iy = ±1 samples as a function of At, 
overlaid with the projection of the likelihood fit result. 
We also performed the CP fit fixing |A| ~ 1, which yields 
sin2/3 = 0.713 ± 0.032 (stat). 

The dominant systematic errors on sin2/3 are due to 
limited knowledge of various background properties, in- 
cuding uncertainties in J/i]j i^-specific backrounds and in 
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FIG. 2: a) Number of r] f = -1 candidates (J/ipK%, ip(2S)K$, 
XciK®, and r} c K%) in the signal region with a B° tag (N B o) 
and with a B° tag (Ng ), and b) the raw asymmetry, (N B o — 
Ng )/(N B o +Ng ), as functions of At. Figures c) and d) are 
the corresponding distributions for the rjf = +1 mode J/ip K®- 
To enhance the signal component, all distributions exclude 
Other-tagged events. The solid (dashed) curves represent the 
fit projections in At for B° (B°) tags. The shaded regions 
represent the estimated background contributions. 

the amounts of peaking backgrounds and their CP asym- 
metries (0.010), to possible differences between the i3fl av 
and Bcp tagging performances (0.009), to the descrip- 
tion of the At resolution functions (0.008), to the knowl- 
edge of the event-by-event beamspot position (0.005). 
The only sizeable sysyematic uncertainties on |A| are 
due to the possible interference between the suppressed 
b — > ucd amplitude with the favored b — > cud amplitude 
for some tag-side B decays [l2j (0.015), and to the CP 
content of the peaking backgrounds (0.006). The total 
systematic error on sin2/3 (|A|) is 0.018 (0.017). The main 
systematic uncertainties on both sin2/? and |A| for the full 
sample, for the seven individual modes, and for the fits 
to the J/ip K° and J/tp K® samples are summarized in the 
Table at 

The large Bcp sample allows a number of consistency 
checks, including separation of the data by decay mode 
and tagging category. The results of those checks, all 
consistent within the errors, are listed in Table Q] We 
observe no statistically significant asymmetry from fits 
to the control samples of non-CP decay modes. 

In summary, we report improved measurements of 
sin2/3 and |A| that supersede our previous results [B|. We 
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measure sin2/3 = 0.714 ± 0.032 (stat) ± 0.018 (syst) and 
| A| = 0.952 ± 0.022 (stat) ± 0.017 (syst), providing an im- 
proved model-independent constraint on the position of 
the apex of the unitarity triangle [l4| . Our measurements 
agree within errors with the published results 15, li[ and 
with the theoretical estimates of the magnitudes of CKM 
matrix elements in the context of the SM 17J . The mea- 
sured value of | A | is consistent with no direct CP violation 
with a significance of 1.72 standard deviations. We re- 
port the first individual measurements of sin2/3 and |A| 
for each of the decay modes within our CP sample, and 
of the J/^K a (K° s + K° L ) sample. 
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